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INTRODUCTION 


I . 

When  a  kinetic  energy  (KE)  projectile  is  shot  against  a  homogeneous  barrier 
(target),  its  performance  is,  in  general,  a  function  that  increases  monotonically 
with  impact  velocity.  The  performance  may  be  taken  as  penetration  depth,  a  limit 
thickness  just  perforated  witn  negligible  residual  mass/velocity  (energy)  behind 
the  barrier,  or  some  other  appropriate  measure  of  effectiveness.  Nearly  all 
physical  penetrator-target  interactions  admit  a  performance  description  in  terms 
of  a  relation  like: 

P=L  •  S(V,pt.  p p.H.  Y,  f,  ...)  ,  (1) 


where  P  is  the  performance,  L  is  a  characteristic  dimension  of  the  penetrator, 
and  s  is  a  characteristic  function  that  depends  on  velocity  V,  the  density  of  the 
target  material  pt,  and  of  the  penetrator  pp,  the  dynamic  strength  of  the  target 
material  H,  and  of  the  penetrator  Y,  the  fineness  ratio  f  (the  length  to  diameter 
ratio  L/D)  and  possibly  other  characteristics  describing  the  interaction.  But, 
S  is  not  explicitly  dependent  on  the  length  or  on  the  diameter  of  the  penetrator. 

The  function  S  increases  with  velocity,  generally  in  a  monotonic  fashion, 
except  in  cases  where  the  character  of  the  penetrator-target  interaction  changes 
in  some  velocity  regime.  For  example,  at  some  velocity,  the  penetrator  may  change 
its  form  rather  suddenly,  say  by  breaking  up,  and  the  S  function  may  dip  in  this 
case  and  then  continue  again  to  increase.  The  formalism  introduced  by  Equation 
1  implies  geometric  scaling  to  be  applicable.  Inter  alia,  this  means  that  the 
same  performance,  P,  can  be  achieved  either  with  a  larger  penetrator  at  a  lower 
velocity  (smaller  S  function  value),  or  a  smaller  penetrator  at  a  higher  velocity 
(larger  S  function  value).  By  applying  additional  constraints,  like  keeping  the 
impact  energy  constant,  it  is  possible  to  show  that  for  most  S  functions,  there 
exists  an  optimum  velocity  where  the  performance  of  the  penetrator,  subject  to 
that  constraint,  is  maximized. 

This  methodology  was  discovered  at  the  BRL  circa  1976  and  has  been  used 
effectively  in  many  studies.  It  has  been  presented  in  "closed"  meetings,  reported 
briefly  on  some  of  its  main  features  (Wright,  19831)  and  its  main  conclusions 
have  been  used  in  commenting  on  the  significance  of  some  recent  hypervelocity 
results  (Silsby,  1984^).  It  was  presented  at  The  1986  Hypervelocity  Impact 
Symposium  held  in  San  Antonio,  TX.  The  proceedings  of  this  symposium  have  been 
published  in  the  International  Journal  of  Impact  Engineering3.  This  report  is  an 
expanded  version  of  the  paper  presented  at  the  1986  symposium. 
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II.  KE  PROJECTILE  PERFORMANCE  SUBJECT  TO  CONSTRAINTS 


One  constraint  of  interest  is  that  of  constant  penetrator  energy.  With  some 
launcher  systems,  like  electromagnetic  guns  (EMG),  it  is  possible,  in  principle, 
to  launch  projectiles  over  a  very  wide  range  of  mass  and  velocity  holding  the 
kinetic  energy  constant,  without  a  significant  impact  on  the  system  size  and 
cost.  The  kinetic  energy  of  a  projectile  is  simply: 


1  nD2 


Lp 


V 2  =  f  D±VL 
p  8/(*Pp 


=  f 


d3v2 
K 3 


(2) 


for  a  cylindrical  projectile  of  diameter  D,  length  L  and  material  density  pp 
moving  with  the  velocity  V.  From  Equation  2,  we  may  calculate  the  diameter  of  a 
constant  energy  projectile  as: 


pl/3 

D=K - - - 

f  l/3y2/3 


(3) 


With  L  =  f  D  and  Equation  1,  we  obtain: 


P=K 


fi/3  £1/3 
V2/3 


S(V,  ...)  . 


(4) 


The  extrema  of  Equation  4  are  found  by  solving  dP/dV=P’=0.  The  equation  for  the 
extrema  of  P,  subject  to  the  constraint  E  =  constant  and  f  =  constant  is: 

V%=VS>  =!  5  for  V=Vopc  .  (5) 


The  solution  of  Equation  5  yields 
the  velocity  (or  velocities)  where 
P'=  0,  and  therefore,  P  assumes  an 
extremum.  For  the  kind  of  S  func¬ 
tions  described  earlier,  this  ex¬ 
tremum  is  a  maximum  of  P  subject 
to  the  constraints  noted  above. 

The  form  of  Equation  5  sug¬ 
gests  a  simple  graphical  means  to 
find  the  optimum  when  the  S  func¬ 
tion  is  not  given  as  an  explicit 
analytic  function  of  V,  but  exists 
only  as  a  graph  of  experimental 
data.  A  line  drawn  from  the  origin  (0,0)  to  the  2/3  value  of  is  parallel 

to  the  tangent  of  the  S-V  curve  at  V  t,  as  illustrated  in  Figure  1. 
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The  application  of  this 
graphical  ta'-'.nique  to  real 
experimental  data  obtained  in 
the  velocity  range  between 
1.2  km/s  and  4.6  km/s  illus¬ 
trates  the  procedure  in  more 
detail.  The  data  of  Silsby 
(19842)  are  summarized  in 
Table  1.  The  S  function  val¬ 
ues  in  the  table  are  calcu¬ 
lated  based  on  solving  Equa¬ 
tion  4  (or  Equation  1)  expli¬ 
citly  for  the  S  function. 

Figure  2  shows  a  plot  of 
these  data  with  a  curve  manu¬ 
ally  fitted  through  the 
points.  Applying  the  above 
technique  to  these  data,  as  shown  in  the  figure,  we  find  with  only  a  few  trials 
an  optimum  velocity  of  about  2.2  km/s.  This  corresponds  to  the  optimum  perf¬ 
ormance  if  the  kinetic  energy  and  the  L/D  are  held  constant.  Note  that  s(vopt> 
...)  is  only  about  1.35,  or  about  85%  of  the  S  value  at  4.5  km/s. 


L/D-23  tungsten  alloy  versus  RHA  data 


VELOCITY  (nv's) 

FIGURE  2  Linear  seals  graphical  fort  bod  for  d»t*rsioln(  1  ha  opt  lsus 

Telocity  If  t  be  kinetic  energy  would  be  held  constant  using 
constant  L/D  tungsten  alloy  eersus  RHA  experiments!  d*i  • 


TABLE  1.  Tungsten  Alloy  Rods  (L/D  =  23,  p=17.3  g/cc)  Versus  RHA 

(8th  International  Symposium  on  Ballistics  -  G.  silsby) 


Mass 

Length 

Diameter 

Velocity 

Depth 

S(V,  ...) 

(q) 

(mm) 

(mm) 

(m/s) 

(mm) 

(P/L) 

100.01 

155.78 

6.87 

1291 

80.0 

0.5135 

98.58 

155.83 

6.82 

1494 

112.0 

0.7187 

98.63 

155.83 

6.83 

1865 

174.5 

1.1198 

46.20 

121.79 

5.28 

2365 

165.2 

1.3564 

47.28 

121.79 

5.35 

2409 

172.3 

1.4147 

96.96 

155.83 

6.76 

2653 

228.5 

1.4663 

95.16 

155.83 

6.70 

2742 

228.0 

1.4631 

43.86 

121.79 

5.15 

2746 

176.3 

1.4476 

97.82 

155.83 

6.80 

3335 

237.9 

1.5267 

96.73 

155.83 

6.70 

3449 

241.5 

1.5498 

46.02 

121.79 

5.27 

3580 

188.6 

1.5486 

94.84 

155.75 

6.69 

4398 

246.9 

1.5852 

99.54 

155.83 

6.86 

4415 

248.1 

1.5921 

45.41 

121.69 

5.24 

4525 

193.7 

1.5918 
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Another  set  of  constraints  that  is  useful  minimizes  the  energy  required  to 
achieve  equal  performance.  By  taking  P  constant  and  f  constant  and  proceeding  in 
the  same  manner  as  before,  we  again  obtain  Equation  5.  The  optimum  velocity 
(representing  the  minimum  energy  to  achieve  the  same  performance  P)  is  found  by 
solving  Equation  5  as  before. 


The  next  case  of  interest  is  that  of  a  projectile  with  the  diameter  and 
kinetic  energy  held  constant.  For  simplicity,  we  will  consider  here  only  the 
special  case  where  S  is  not  dependent  on  f.  Calculating  f  from  the  projectile 
energy: 


K'E 

d3v2 


(6) 


and  using: 


P  =f  D  S{V, 


.),xLLflsiv, 


(’) 


and  setting  dP/3v  =  P'  =0,  we  obtain  the  extrema  condition  to  be: 

V  S'  =  2  S  for  V  =  Vopt  .  (8) 

A  graphical  solution 
analogous  to  that  outlined 
for  Equation  5  can  be  devel¬ 
oped  (as  illustrated  in  Fig¬ 
ure  3)  but  not  all  S  curves 
yield  a  solution.  A  unique 
solution  requires  an  S-V 
curve  with  a  relatively  steep 
slope. 

An  even  simpler  geometric 
construction  is  possible  if 
the  data  are  plotted  using 
logarithmic  scales,  i.e, 

£n(P)  vs.  £n(V).  Since  the 
log  transformation  is  mono¬ 
tonic,  the  extrema  of  P  vs.  V  occur  at  the  same  values  as  those  of  £n(P)  vs. 
£n(V).  Following  the  same  procedures  as  used  above,  we  find  that  maximum 
performance  at  constant  energy  and  constant  L/D  occurs  when  the  d (In  P)/3(ln  V) 
=2/3  and  the  maximum  for  constant  energy  and  constant  diameter  occurs  when  3(£n 
P)/d(£n  V)  =  2.  In  all  cases,  the  underlying  assumption  is  that  geometric  scaling 
applies. 


3 
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VELOCITY 

Method  for  determining  the  opt  Ima  velocity  (or  •  given  9 
function  when  kinetic  energy  end  diameter  ere  be  Id  constant. 
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An  example  of  the  graph¬ 
ical  method  of  finding  maxi¬ 
mum  performance  for  constant 
energy  and  constant  L/D  using 
loqarithmic  scales  is  shown 
in  Figure  4.  The  data  are  the 
same  as  that  of  Figure  2  and 
the  curve  is  manually  fitted 
through  the  points.  A  refer¬ 
ence  line  which  has  a  slope 
of  2/3  is  shown  and  a  line 
parallel  to  it  but  tangent  to 
the  data  curve  yields  the 
point  of  maximum  performance 
if  the  kinetic  energy  had 
been  held  constant.  Again, 
the  velocity  corresponding  to 
this  point  is  about  2.2  km/s. 

The  methodology  outlined  above  is  most  general  in  nature.  The  S  curves  may 
be  analytic  forms,  may  be  given  in  graphical  form  based  on  experimental  data  or 
may  be  generated  by  numerical  methods. 


FIGURE  4  Logarithmic  acale  graphical  method  for  determining 
the  opt  ImuiD  velocity  If  the  kinetic  energy  would  be 
held  constant  using  constant  L/D  tungsten  alloy  versus 
RHA  experimental  data 


III.  PENETRATION  OF  HOMOGENEOUS  TARGETS 

Next,  we  shall  provide  some  examples  using  explicit  analytic  formulations 
that  are  suitable  throughout  most  of  the  hyper-velocity  regime.  For  simplicity, 
we  will  restrict  the  discussion  to  the  most  simple  case  -  penetration  of 
homogeneous  targets.  First,  we  deal  with  the  short  penetrator  of  constant  mass. 
Vitman  and  Zlatin  (1962)4  and  Belyakov,  Vitman  and  Zlatin  (19625,  19636)  have 
shown  that  the  resistance  to  penetration  is  given  by: 

H  +  kg  pt  V2  ,  (9) 

where  H  is  a  dynamic  material  property  which  is  only  weakly  dependent  on 
velocity  They  showed  that  the  resulting  simple  relation: 

P= L-- ^  hix*2k0zl)  where  z£=  9 ^  (1°) 

2k0  pc'  2 H 

and  Vq  is  the  striking  velocity,  describes  the  penetration  behavior  of  different 
penetrator  and  target  materials  rather  well  to  very  high  velocities.  In  the 
following,  we  take  2kg  =  1  and  add  to  Equation  9  a  term  that  accounts  for  the 
effect  of  the  target  material  moving  in  front  of  the  penetrator  in  a  non-steady 
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flow  (Goodier,  19647;  Hanagud  and  Robs,  19718;  Byrnside,  Torvik  and  Swift, 
19729;  Johnson,  197210;  Tate,  197911).  In  passing,  we  note  that  the  target 
resistance  term  is  directly  related  to  the  work  required  to  make  a  unit  volume 
cavity  in  the  target  (Johnson,  197210;  Hill,  198012;  Hopkins,  I96013).  For 
plastically  deforming  materials,  a  good  estimate  for  it  can  be  obtained  from 
bi-linear  approximations  to  quasi-static  stress-strain  data  (Hill,  198012; 
Hopkins,  I96013)  . 


The  resulting  relation  for  a  non-deforming  constant  mass  penetrator  is  given 
by  (Goodier,  19647) : 


jgo  .Z) 
L 


(11) 


where  the  quantity  ptD/(3ppL)  accounts  for  the  additional  accelerated  target  mass 
in  front  of  the  penetrator  (the  factor  1/3  implies  a  hemispherical  cap  of 
affected  target  material  in  front  of  the  penetrator).  Note  that  Z2  =  ptV2/(2H). 
(See  Appendix  A  for  derivation  of  Equation  11.) 


Setting  S(ZQ,Z)  equal  to  the  right  side  of  Equation  11  (see  Equation  1), 
setting  Z  =  0,  taking  the  derivative  of  S(ZQ,0)  with  respect  to  ZQ  and  making 
appropriate  substitutions  into  Equation  5  (replacing  V  by  Z  since  Z  is  directly 
proportional  to  V),  the  optimum  ZQ  at  constant  kinetic  energy  is  found  by  solving 

(l+Zo)  «n(l+Zo)-3Zo  =0  (12) 


which  has  the  approximate  solution  of  Zn  =  3.975. 


For  a  target  material  of 
armor  steel,  the  value  of 
2H/pt  is  about  (1000  m/s)2, 
so  that  optimum  penetration 
performance  occurs  at  an  im¬ 
pact  velocity  near  4000  m/s. 
The  curve  in  the  vicinity  of 
the  optimum  velocity  is  near¬ 
ly  flat.  Therefore,  near  op¬ 
timum  penetration  for  con¬ 
stant  impact  energy  can  be 
achieved  over  a  wide  range  of 
velocities.  This  is  illustra¬ 
ted  in  Figure  5  for  a  one 
megajoule  L/D=2  iron  penetra¬ 
tor  impacting  armor  steel. 


H 

2 


S  2  I  0.5 


MASS 

0.1 


(kg) 


t  -  I  Ul  H  :  S  6  GP.  r  =  z 

Pp  =Pt=  7«50 


:  5* 


0  Z  4  6  B  10 

IMPACT  VELOCITY  (km/s) 

riCURE  5  Penetration  performance  of  an  iron  projectile  striking 
a  ateel  target  with  one  mega-Joule  of  kinetic  energy 
showing  optimum  performance  and  5«  below  optimum 


12 


Of  much  greater  interest  in  practice,  are  deforming  or  eroding  penetrators. 
There  have  been  many  papers  published  on  this  subject  beginning  with  the 
classical  works  of  Alekseevskii  and  Tate  in  the  1960's  (Alekseevskii,  1966^; 
Tate,  196715, 196916)  with  more  refinements  which  remove  some  of  the  earlier 
limitations  being  added  even  today  {Tate,  198617) .  In  the  following,  we  shall 
present  formulations  (using  a  consistent  notation)  which  we  have  found  useful. 

Following  the  basic  ideas  of  Christman  and  Gehring  (19661®),  the  penetration 
of  a  long  deforming/eroding  rod  penetrator  may  be  separated  according  to  the 
important  identifiable  phases  or  stages  of  the  penetration  process.  For  example, 

stage  1  stage  2  stage  3 
P(Z0.Z)=m1D  SD1(Z0 ,  Z^+^L—m^D}  S^Zy,  Z^+in^D  S^jZ^tZ^  (13) 


where  ZQ>Z1>Z2>Z>0.  The  first  term  describes  the  transient  where  m1  diameters  of 
the  rod  are  "used"  up.  The  second  term  is  the  quasi-steady  state  phase  which 
erodes  the  length  L  -  mjD  of  the  rod.  The  final  stage  is  the  penetration  by  a  rrijD 
constant  mass  penetrator.  It  is  not  obvious  that  the  equation  described  above  can 
always  be  constructed  in  a  consistent  manner.  However,  we  have  found  that  the 
following  simpler  form  gives  satisfactory  results: 

P(Z0,Z)  =  (L  -  mD)  Sj(Z0,Z)  +  mD  SD(ZQ,Z)  .  (14) 

This  may  be  considered  an  approximation,  in  the  spirit  of  Christman  and  Gehring, 
which  holds  over  the  whole  velocity  range,  including  low  velocities  where  their 
original  formulation  fails.  The  diameter  function  SD,  is  based  on  Equation  11: 


SD(Zo  > 


(15) 


Here,  m  is  considered  an  adjustable  parameter  with  typical  values  between  0.5 
and  2.  In  this  formulation,  SD  attempts  to  describe  both  stage  1  and  stage  3  of 
the  penetration  process.  When  used  by  itself  (non-eroding  long  rod  penetrator), 
m  is  set  equal  to  L/D.  For  the  S  function  denoted  as  Sj(Z0,Z),  we  will  present 
several  functions  appropriate  to  special  cases  of  interest.  (More  details 
concerning  these  S  functions  are  in  Appendix  B.) 


The  first  function,  S1f  is  that  for  a  penetrator  with  no  strength  (  Y  -*  0, 
using  the  Tate/Alekseevskii  notation)  representing,  for  example,  a  non-stretching 
jet  (Rostoker,  195319) .  Since  such  a  penetrator  does  not  slow  down  during 
penetration,  S  is  a  function  of  ZQ  only,  and  the  process  terminates  when  the 
length  L  -  m  D  is  used  up. 
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(16) 


«i(Zo)-' 


3V/y02+i-n2 


n(l+202) 


for  Zc <  Z0  s  » 


and  Sx(Z0)=0  for  0<Z0<Zc=M/  where  a2  =  pt/pP  (Tate/Alekseevskii  notation)  and 
Zc  is  the  normalized  cut-off  velocity  below  which  no  penetration  occurs.  See 
Appendix  B  for  some  background  on  the  derivation. 


The  second  S  function  for  the  length  term  is  Tate's  special  case  for  Y  =  H, 
applicable  to  the  entire  velocity  range  from  0  to  infinity.  (See  Appendix  B  for 
some  background  on  the  derivation. ) 


S2{Zo'Z)  =  ± 


(17) 


If  we  subject  S2(Z0,0)  to  the  optimization  procedure  for  constant  energy  and  L/D 
,  i.e.,  solve  Equation  5,  we  obtain  tne  transcendental  equation: 

8n(l+3y)-y=0  ,  (18) 


where 


y= 


-•opt 


p(l  +  )i) 


which  has  the  solution  y  =  1.9038.  With  /i  <  1  (high  density 


0  w 

penetrator  impacting  lower  density  targets),  the  optimum  velocity  Vopt  =  Zopt  ^  ~ 
can  be  considerably  lower  than  that  for  the  case  of  the  constant  mass  penetrator 


discussed  earlier. ‘ 


Tate  also  gave  a  solution  for  Y  equal  to  fractional  values  of  H,  but  only 
for  the  restricted  case  of  p=l  (target  and  penetrator  densities  equal).  Using  an 
approximate  solution  of  the  Tate/Alekseevskii  "Bernoulli  equation",  but  with  a 
better  approximation  than  that  of  Alekseevskii ' s,  that  preserves  the  cut-off 
velocity,  we  developed  a  closed  form  solution  for  the  case  Y  =  H/(2+p),  which  is: 


Hz< 


’  n  a0 


(19) 


for  ac<a0£  °°  and  S3(Z0,0)=0  for  0<a0<ac  where  a0  =  Z2|  —  |  and  a c  =  fj.  See 

Appendix  B  for  the  derivation. 


More  functions  of  this  kind  can  be  given,  but  only  for  Y  values  that  are 
smaller  than  that  for  S3.  Since  the  S3  and  S3  curves  are  not  greatly  separated, 
these  new  functions  are  not  of  great  practical  value.  What  is  needed,  however, 
is  a  family  of  functions  that  till  the  large  gap  between  S2  and  SD.  Functions  of 
this  kind  can  be  generated  by  using  what  must  be  considered  a  radical,  if  not 


14 


bizarre,  approximation  to  the  variable  mass  penetrator  problem. 

In  experiments  with  high  density  long  rod  penetrators  against  finite 
thickness  plates,  it  is  noted  that  the  length  of  rod  lost  and  eroded  is  pro¬ 
portional  to  the  target  thickness  perforated.  If  systematic  variations  with 
velocity  are  ignored,  we  can  start  with: 

L*  =  L*  -  rx  =  L0  -  mD  -  rx  (20) 

where  Lx  is  the  instantaneous  rod  length,  L0  is  the  initial  length,  D  is 
penetrator  diameter,  x  is  the  penetration  depth  and  m  and  r  are  dimensionless 
constants.  Differentiating  Equation  20  and  noting  that  dL/dt  =  -(v-u)  where 
(using  Tate's  notation)  v  is  the  instantaneous  velocity  of  the  rear  of  the 
penetrator  and  u  =  dx/dt  is  the  penetrator-target  interface  velocity,  we  find 
that  r  =  (v-u)/u.  (This  will  be  referred  to  as  the  constant  v/u  approximation.) 


Next,  the  kinematics  of  the  eroding  rod  are  approximated  by  using  rite's 
"Bernoulli  equation"  and  noting  that  Y  =  ~PpL  dv/dt  (see  Appendix  B) : 

-PpL^  +  ±pp{v-u)**H+±p  tu*  .  (21) 


This  is  now  interpreted  as  an  equation  of  motion  and,  as  with  the  SD  function, 
supplemented  on  the  right-hand  side  with  the  added  target  mass  term: 


■ppidF+ip^<v'u)2s'w+^PtU2+^Dpt'^ 


dt 


(22) 


With  the  earlier  constant  v/u  approximation,  we  have  v  =  u(l+r),  (v-u)2  =  (ru)2 
and  "5^ =  •  Therefore,  Equation  22  can  be  integrated  after  dividing 
through  by  pp  and  substituting  for  L  (using  Equation  20),  (v-u)2  and  du/dt  and 
then  separating  variables.  By  defining: 


c2  =  -i* 


H2(l+r2) 


0<  c2  < 1 


(23) 


fc.iiilil  >0, 


ji2-r2 


(24) 


c-  P_tg  -  U2 

3 P pL„  3p p{L0-mD)  3  (f-m) 


and 


d= 


e 

1+r  ' 


the  result  is: 


(25) 


(26) 
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The  primary  parameter  is  r  which  is  the  proportionality  factor  between  rod  mass 
loss  and  penetration.  It  can  be  shown  that  as  r  -*■  p,  S4  -*  S2.  Also,  as  r  -*  0,  the 
differential  equation  (Equation  22)  approaches  that  for  SD.  Therefore,  the  range 
0<  r  <p  generates  curves  between  SD  and  S2.  The  physical  relevance  of  Equation 
27  can  not  be  assessed  at  this  time,  so  it  should  not  be  taken  too  seriously.  The 
constant  v/u  approximation  which  is  the  basis  for  deriving  S4  is  certainly  open 
to  severe  criticisms.  Nevertheless,  the  collection  of  S  functions  given  above, 
together  with  Equation  14,  allows  for  some  interesting  parametric  studies. 

The  character  of  these  functions  is  shown  in  Figure  6.  The  functions  Si,  S2, 
S3  and  S4  (with  r«p)  tend  to  the  high  velocity  limit  of  1/p  while  S4  with  r=0 
approaches  the  SD  function.  The  S„  function  tends  toward  infinity  in  a  loga¬ 
rithmic  fashion.  See  Appendix  C  for  some  additional  examples  of  parametric 
variation. 
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IV.  APPLICATIONS 


An  example  of  the  appli¬ 
cation  of  these  analytical 
forms  is  shown  in  Figure  7 . 

This  is  a  plot  of  the  final 
penetration  P  (Equation  14 
with  Sj  being  replaced  by  S1 
through  for  each  curve 

respectively)  as  a  function 
of  the  impact  velocity  for 
the  case  where  E=1  MJ,  f=20, 
p=0.7,  r=p/2,  m=l  and  pt  = 

7850  kg/m*  (armor  steel).  The 
performance  (note  that  the 
ordinate  is  P  and  not  P/L) 
maxima  (the  peak  of  each 
curve)  occur  at  P  =  0.30, 

0.41,  0.34,  and  0.64  meters 
with  the  optimum  velocities 
of  2.0,  1.5,  1.72  and  1.52  km/s  for  the  functions  S2,  S3  and  S4  respectively. 
The  non-dimensional  Z  values  are  2.09,  1.57,  1.79  and  1.59  respectively.  The 
penetration  at  the  optimum  velocity  increases  with  penetrator  strength  Y  since 
the  functions  S1<S3<S2<S^  are  associated  with  the  penetrator  strengths  0  = 
Y^<Y3<Y2<Y^,  and  the  corresponding  cut-off  velocities  decrease.  If  the  dynamic 
strength  of  the  target  material,  H,  is  increased,  the  performance  maxima 
(indicated  in  Figure  7  by  a  tic  mark  at  the  peak  of  each  curve)  decrease  in  value 
and  shift  to  higher  velocities. 

The  curves  in  Figure  7  are  plots  of  Equation  14  with  E  =  1  MJ.  If  the 
ordinate  of  the  graph  is  re-scaled  and  labeled  as  P/(E1/3)  in  units  of  m/(MJ)1/3, 
Figure  7  may  be  used  for  other  constant  energy  levels  as  well.  For  example,  with 
E  =  8  MJ  or  27  MJ,  the  P-values  of  Figure  7  would  be  2  times  or  3  times  higher 
respectively. 

The  experimental  data  from  Table  1  (Section  II)  is  for  a  tungsten  alloy  (WA) 
penetrator  fired  against  RHA.  The  dynamic  strength  for  WA  is  between  1  and  2  GPa 
and  the  density  is  17300  kg/m3.  The  dynamic  strength  of  RHA  is  between  4  and  6 
GPa  and  its  density  is  7850  kg/m3.  For  these  densities,  p  =  y/p  t/pp  =  /0 . 4537  =0.6736  . 
The  criterion  for  developing  the  S3  function  is  that  Y  =  H/(2  +  p) .  Substituting 
the  values  of  4  and  6  GPa  for  H  and  the  value  for  p  yields  Y  =  1.5  and  2.2  GPa 
respectively.  Applying  Equation  14  with  S3  for  Sjf  taking  m  =  1,  kQ  =  1,  k{  =  kpf 


MASS  (kg) 


FIGURE  7  Two  phase  penetration  performance  (  f  -  m)  Sj(Zo.O)  * 
mD  Sd(Zq,0).  The  mark  at  the  peak  of  each  curve 
shows  the  velocity  for  optimum  penetration. 
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(see  Appendix  B)  and 
f=23,  the  curves  of 
Figure  8  are  pro¬ 
duced  (taking  H  =  4, 

6  or  8  GPa).  The 
data  points  are 
those  from  Table  1. 

Above  1800  m/s,  the 
data  lie  between  the 
curves  for  H  =  4  and 
H  =  6  GPa.  The  two 
points  below  1800 
m/s  lie  between  the 
H  =  6  and  H  =  8  GPa 
curves.  These  two 
shots  were  fired  at 
the  BRL  while  the 
rest  were  fired  at 
the  hypervelocity 
range  facility  at  Tullahoma,  TN.  Whether  the  inconsistency  is  due  to  experimental 
differences  or  that  the  criterion  for  the  Sj  function  is  not  applicable  is  not 
known  at  this  time.  It  seems  certain  that  a  value  of  8  GPa  for  the  dynamic 
strength  of  RHA  is  much  too  high. 


There  are  other 
way 8  of  getting  a 
curve  to  fit  the  data. 
For  example,  shifting 
the  S3  function  to  the 
right  by  390  m/s,  tak¬ 
ing  m  =  1.404,  and  H  = 
4  GPa  results  in  the 
curve  of  Figure  9.  It 
is  not  clear  that 
shifting  the  S3  func¬ 
tion  is  justified. 
What  is  needed  is  a 
set  of  data  for  one  or 
more  penetrators  for 
striking  velocities 
from  about  200  m/s  out 
to  4500  m/s  at  inter¬ 
vals  less  than  500  m/s 


FIGURE  9  Fitting  the  data  of  Table  1  using  the  two  phase 

penetration  performance  of  Sp(Zo.O)  and  a  shifted 
S3(Zq.0)  with  H  =  4  GPa 


FIGURE  8  Bracketing  the  data  of  Table  1  using  the  two  phase 
penetration  performance  of  Sd(Zo,0)  and  S3(Zo,0) 
with  H  =  4,  6  and  8  GPa. 
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where  both  Y  and  H  are  known  to  a  good  degree  of  certainty. 


In  applying  these  equations  to  a  particular  situation,  one  should  keep  in 
mind  that  increasing  the  value  of  H  will  shift  the  curve  to  the  right.  Increasing 
the  value  of  pt  shifts  the  curve  down.  Increasing  the  value  of  Y  will  shift  the 
curve  to  the  left.  Increasing  pp  shifts  the  curve  up.  In  addition,  one  should  be 
cognizant  of  the  following  which  are  based  on  experimental  experience  with  metal 
impacting  metal  interactions. 


(1)  If  Y  <  H  the  penetrator  will  erode  during  penetration 

and  for  v(s7c  =J  ~  —  (see  Appendix  B,  equation  B-6)  the  depth  of 


N  kpPP 

penetration  is  only  a  few  milli-meters. 


(2)  If  Y  >  H,  the  penetrator  remains  rigid  -  doesn’t 

deform  unless  v>.  where  V  is  the  instantaneous  velocity. 

N  ktPt 

(3)  For  a  given  set  of  conditions  which  include  normal 
impact,  a  penetrator  will  defeat  a  finite  thickness 
target  with  a  substantial  residual  velocity  when  the 
target  thickness  equals  the  depth  of  penetration  into 
the  semi-inf inte  target. 


V.  CONCLUDING  REMARKS 


These  examples  illustrate  some  of  the  capabilities  of  the  methodology 
presented  in  this  paper.  Another  obvious  application  is  the  first  order  analysis 
of  multiple  penetrators,  that  is,  the  performance  of  one  long  penetrator  may  be 
compared  to  N  shorter  penetrators  of  equal  material  length  and  mass  (segmented 
penetrators).  This  last  case  can  not  consider  the  effects  of  separation  between 
penetrator  segments,  nor  can  it  account  for  penetrator  mass  remaining  in  the 
crater  from  previous  impacts  in  a  successive  series  of  impacts.  Never-the-less, 
it  provides  for  the  initial  assessment  of  the  efficacy  of  segmented  penetrators, 
especially  in  the  higher  velocity  regime. 

Finally,  we  will  attempt  to  answer  one  obvious  question.  In  this  day  and  age 
of  "super  computers"  and  refined  computational  techniques,  why  bothei  with 
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approximate  closed  form  analytic  expressions?  In  the  exploratory  stages  of 
research  and  development,  there  are  many  occasions  where  a  large  number  of 
parametric  investigations  need  to  be  made.  The  requirement  is  not  necessarily  the 
non-plus-ultra  in  accuracy,  but  a  consistent  set  of  equations  that  reflects  the 
major  trends  reasonably  well  and  does  not  misrepresent  the  important  physical 
interdependencies.  For  this  kind  of  use,  the  formulations  presented  in  this  paper 
prove  to  be  superior.  At  their  best,  they  provide  reasonable  first  estimates 
requiring  a  minimal  amount  of  material  parameters  and  of  computational  resources. 
At  their  worst,  one  can  delegate  them  into  the  realm  of  semi-empirical  fitting 
formulas,  but  with  parameters  that  have  some  direct  physical  relation  to  the 
penetrator-target  interaction  process. 
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APPENDIX  A 

SHORT  NON-ERODING  (CONSTANT  MASS)  PENETRATORS 
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Ideally,  the  relations  listed  below  apply  (in  a  strict  sense)  only  to  a 
rigid  body  penetrator,  i.e.  a  penetrator  that  does  not  change  its  shape  or  its 
mass  during  the  penetration  phase.  With  minor  misgivings,  the  equations  are 
applicable  to  a  deforming  penetrator  provided  it  is  sufficiently  short  so  that 
the  major  portion  of  the  mass  which  has  expanded  radially  retains  a  forward  axial 
velocity.  Taking  Equation  9  of  Section  III: 


H  +  k^V2  (A-l) 

as  the  target  resistance,  letting  =  1/2  and  adding  a  term  which  accounts  for 
non-steady  flow,  we  find  the  following  expression  for  the  pressure  at  the 
penetrator-target  interface: 

Pt  =  W+^P  t^2+4p  •  (A-2) 

2  3  at 

This  pressure  is  distributed  over  the  presented  area,  A  =  trD2/4,  of  the 
penetrator  which  is  assumed  to  be  constant  or,  at  most,  slowly  varying.  The  mass, 
M  =  A  L  pp,  is  taken  to  be  constant.  To  a  first  approximation,  the  equation  of 
motion  for  the  penetrator-target  system  follows  from  balancing  the  axial  forces: 

Hjg.ALPpVjg~Ap,  (A-31 


which  can  be  integrated  in  a  close',  form  solution  to  produce  Equation  11  of 
Section  III.  This  equation  is: 


^Zo.Z)  =  L^  1* 


(A-4) 


where  Z2  =  ptv2/(2H)  and  P(ZQ,Z)  is  the  penetration  depth  when  the  instantaneous 
velocity  is  V.  ZQ  is  Z  with  V  equal  to  the  striking  velocity  V0.  The  accompanying 
figures  illustrate  the  penetrator-target  intaraction  described  in  this  Appendix. 


Equation  A-4  is  not  applicable  to  situations  where  the  shape  factor  kQ  is 
not  constant  during  the  penetration  interaction.  The  value  used  here,  1/2,  is 
applicable  to  a  blunt  nosed  penetrator. 
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APPENDIX  B 

THE  ALEKSEEVSKII-TATE  MODEL  FOR  ERODING  LONG  ROD  PENETRATORS 
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B.l  Background 


The  one-dimensional  eroding  Long  Rod  Penetration  Model  as  formulated 
independently  by  Alekseevskii  and  Tate  is  based  on  the  following  simplifying 
assumptions: 

1.  The  dynamic  behaviour  of  the  homogeneous  penetrator  is  that  of  a  rigid 
body  up  to  some  pressure  Y.  Likewise,  the  homogeneous  target  acts  as 

a  rigid  body  up  to  some  pressure  H.  The  material  flows  like  a  fluid 
for  pressures  above  Y  for  the  penetrator  or  H  for  the  target. 

2.  After  an  initial  impact  transient,  the  penetrator-target  interface  moves 
with  a  velocity  U.  If  the  penetrator  erodes  (becomes  shorter  than 

the  initial  length  LQ  while  penetrating  the  target)  the  tail  end  of 
the  rigid  portion  of  the  penetrator  has  a  velocity  V  while  the  front 
end  has  a  velocity  V  -  U  relative  to  the  interface. 

3.  The  retardation  of  the  rigid  portion  of  the  long  rod  penetrator  (LRP)  of 
length  L  is  determined  only  by  the  pressure  Y  -  the  maximum  stress 
level  the  still  rigid  rod  material  can  sustain  before  it  begins  to 

act  as  a  fluid. 

The  statements  above  are  summarized  by  the  following  coupled  differential 
equations  (refer  to  Figure  B-l  with  respect  to  Equations  B-l  and  B-2): 


(B-l) 


~  =  -(V-U)  ,  (B-2) 


ApLpp^u  (B'3) 

In  these  equations,  P  is  the  depth  of 
penetration,  L  is  the  instantaneous 
length,  Ap  is  the  penetrator  cross- 
sectional  area,  and  pp  is  the  penetra- 

Up  to  this  point,  the  Alekseevskii  and  the  Tate  formulations  are  identical. 
However,  there  are  only  three  equations  with  four  unknowns  (P,  L,  U  and  V).  So 


| GIVEN  A  TIME  INCREMENT  At  \ 

DISTANCE  TRAVELLED  WITH  NO  TARGET 
-  V  At  - 


LENGTH  ERODED 
- (V-U)  At  — 


-  aL 


DISTANCE  TRAVELLED  IN  TARGET 

-  U  At  - 

-  aP  - 


THEREFORE:  s  -  (v-u)  ®„<t  =  u 

At  At 


FIGURE  B-l  Basis  for  dL/dt  and  dP/dt  equations. 


tor  density. 


B  -  3 


an  additional  relation  is  needed  to  solve 
the  system  of  equations. 


Alekseevskii  choses  this  relation  to  be 
the  force  (or  pressure)  balance  across  the 
penetrator-target  interface: 

Ap  (y+Jcpp  p{V-U)2)  =  Ap  (H+k'P  CU 2)  ( B-4 ) 


>  = u  where  kp  and  kt  are  shape  factors  which 

•  ^ ^  MFT|| — *  characterize  the  deformed  regions  in  the 

penetrator  and  target  materials  near  their 
r,CU“B'2  respective  rigid-fluid  boundaries.  In  gener- 

_ al,  both  of  these  factors  are  approximately 

1/2. 

Tate,  on  the  other  hand,  uses  a  modified  "Bernoulli  Equation"  for  the 
stagnation  pressure  of  two  colliding  fluid  streams: 


Y+±pp(V-U)2  =  ±ptU2+R  . 


(B-5) 


where  Y  and  R  have  the  same  meaning  as  Alekseevskii ' s  Y  and  H.  When  kp=kt=l/2. 
Equations  B-4  and  B-5  are  identical.  However,  their  physical  justification  is 
different. 


target  surface 

RICIO  PORTION  OF 
PENETRATOR 
<  •  f  I  »D|  »k  L> 
TRAVELLING  AT 
VELOCITY  V 


penetrator  target  interface 
MOVING  uith  velocity  u 


AT  INTERFACE 


(1H0  of  ar*ou«>  *  i 
■ - : - : - r~  »  ■» 

♦  Y  ♦  ♦  *  1 


PRESSURE  AT 
INTERFACE 
<CHO  OF  —out) 


PENETRATOR  MATERIAL 
EFFLUENT 


Y.-PL^ 
P  6  ♦ 


AXIAL  PRESSURE  BALANCE  AT  PENETR A TOR -TARGET  INTERFACE 
< AVERAGED  OVER  PENETRATOR  CROSS-SECTIONAL  AREA) 

y  *  kp  Pp<  V  -  U)S  ■  H  ♦  k.f.U2 


FIGURE  B  ~3  Eroding  p»n»  *ro  tor  i  n  tor o<  t  i  ng  with  forget 
under  Steady  ptni  tro  t  ion  conditions. 


Equation  B-4  (or  Equation  B-5)  delimits  the  special  case  of  rigid  target  or 
of  rigid  penetrator  behaviour: 

a.  Rigid  target  defined  by  U  =  0,  H  2  Y  and  V2 <:  V%  =  Y  .  (B-6) 

*pPp 

b.  Rigid  penetrator  defined  by  (V  -  U)  =  0,  H  <  Y  and  V2  =  U2  <.  .  (B-7) 

^tP  t 


The  regime  of  interest  is  that  of  the  eroding  penetrator  and  the  target 
being  penetrated,  namely,  (V  -  U)  >0  and  U  >  0.  This  corresponds  to  the  velocity 
range  V£  s  V  <  <*>  or,  equivalently,  the  parameter  range  1  a  a  >  0. 


Setting: 


g.  H-Y  _V*c 
kpPpV2  V2 


and 


t 

^pPp 


(B-8,9) 


the  exact  solution  for  U  of  the  quadratic  equation  (Equation  B-4)  taking  the 
negative  sign  of  the  radical  is: 


U=  V 


1-fl _ 

l  +  Vn2+«  (1-p2) 


(B-10) 


(Note  that  a  is  not  constant,  i.e.,  it  iB  an  inverse  function  of  V.)  To 
facilitate  integrating  Equations  B-l  through  B-3,  it  is  necessary  to  use  an 
approximate  relation  for  U.  Alekseevskii  uses  the  approximate  solution: 


U _ 1_ _ a_ 

V  l  +  n  2)1 


(B-ll) 


which  diverges  substantially  from  the  exact  solution  near  the  cut-off  velocity 
Vc  unless  n*l.  The  cutoff  velocity  occurs  when  a  =  1  and  where  penetration 
ceases.  At  higher  velocities,  it  is  quite  satisfactory.  However,  a  better 
approximation  is 


U  _  l-«  V-U  _  a+\l 

V  l  +  ^  "*  V  1  +  n 


(B-12 ) 


which  yields  the  correct  limit  of  U  =  0  when  V  =  Vc  (that  is,  a  =  1)  and  which 
differs  by  only  a  few  percent  from  the  more  complicated  exact  solution  (Equation 
B-10)  over  the  velocity  range  of  interest:  V  >  Vc. 


B.2  Erosion  (change  in  length) 

Combining  Equations  B-2  and  B-3,  the  differential  equation  for  the  change 
in  length  can  be  expressed  as: 


B-5 


(B-13) 


dL=^JL  (V-U)  dV 


Substituting  the  approximation  (Equation  B-12)  for  (V  -  U) s 

|2  VdV  . 


Setting: 


dL~^-£L 


2  Y 


JilfLl 

i+»* 


t)  =  —  whereby  2VdV=V%dr\  ; 


«o  = 


H-Y 

kpPpVt 


Jl 

Vo 


a  = 


n  = 


Ppvj 


2Y(l+n) 


and 


2k  n-  H-Y 

-  m* 1» 


- 


l+2ic_n  (1+n) 


we  obtain 


dL 

L 


Integrating  this  equation  with  the  initial  conditions  (L  =  LQ,  V  =  VQ,  q  = 
the  condition  (L(l,r;),  V,  q)  at  a  later  stage,  the  result  is: 

■£<(1/1))  _„n  ,-a0(l-n) 

- - — - i|  e 

"0 


( 


This  equation  is  valid  for  all  velocities  between  Vc  and  VQ,  that  is,  V( 
VQ.  At  the  cut-off  velocity  V  =  Vc  (when  U  =  0),  the  cut-off  value  for  r 
=  nfi/aQ  and  the  length  which  remains  (not  eroded)  iB: 

L(l,r\c) 


(/Jl»-a0) 


Equations  B-21  and  B-22  are  valid  for  all  values  of  n. 


B-14 ) 

B-15 ) 

B-16) 

B-17 ) 

B-18 ) 

B-19 ) 

B-20 ) 
1)  and 

B-21) 

.  <  V  < 
7  is  lc 

B-22) 
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B.3  Penetration 


Combining  Equations  B-l  and  B-2  yields  the  differential  equation  for  the 
change  in  penetration  depth  as  a  function  of  change  of  penetrator  lengths 

dP=  dL  «  ZSIZULL  dL  .  (B-23 ) 

V-U  \x+a 


As  a  function  of  change  of  velocity  (using  Equation  B-14),  the  change  in 
penetration  depth  can  be  expressed  as  follows: 


dP=  - 


2  y 


l-q 
1  +  M 


2  VdV  . 


(B-24) 


Using  Equations  B-15  ,  B-17  and  B-16,  this  becomes: 

dP~LT  (l-T)dn  ' 


(B-2  5) 


Substituting  L(1,j7)  of  Equation  B-21  for  L  and  using  Equation  B-18,  Equation  B-25 
can  be  re-written  as: 

-f-iH-t  =  -e-a°  f  |-£°_  dT1  .  (  B-26 ) 

A>  J  v  ^  n/ 

By  adding  and  subtracting  ntln_1ea|1''  and  multiplying  both  sides  by  jj,  Equation 
B-26  can  be  partially  integrated  to  yield: 

=Constant-x\  /Jea°<n"1)  +^(i  +  fi)  e~a°  fr)11'1  e*0*1  dr)  .  (B-27) 

L0  J 


The  constant  is  evaluated  by  taking  the  initial  conditions  to  be  (V=V0,  q= 1, 
L=L0,  P=0)  and  integrating  to  (V,  q,  L,  P(l,r;);  q>qc) .  This  results  in: 


H  P(-1>nI  =  i-n  «  e-aod-n)  +  n  (1  +  ^  p111 ' 

Ln  J  i 


n^e^cfn 


(B-28) 


where  qc  =  np/a0  and  P(q<qc)  =  0.  Note  that  the  second  term  to  the  right  of  the 
equal  symbol  is  the  Bame  as  Equation  B-21. 

Closed  form  solutions  for  integer  values  of  n  (n  =  0,1,2,...)  can  be  found. 
The  expressions  for  P(l,rj)  become  progressively  more  complex  as  the  value  of  n 
increases.  The  following  segment  evaluates  the  closed  form  solutions  for  n=0,l 
and  2 . 
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The  case  n  =  ®  corresponds  to  Y  =  0,  an  LRP  without  any  strength,  which  is 
represented  by  a  "fluid  rod".  Equations  B-l  through  B-4  still  apply  but  a 
different  approach  to  the  solution  is  required.  Since  Y  =  0,  there  is  no 
retardation  of  the  portion  of  the  rod  which  has  not  impinged  on  the  target.  So 
the  velocity  of  the  rod  remains  constant  at  V0.  The  interface  velocity  U  is  still 
found  from  Equation  B-12  but  it  also  remains  at  a  constant  value  U0.  However,  its 
value  is  non-linearly  dependent  on  V0.  Integrating  Equations  B-l  and  B-2  with  the 
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(P-34) 


initial  conditions  at  t  =  0,  P(t)  =  0  and  L(t)  =  L0: 

Pit)  =  U0t  . 


Lit)  =L0-  iV0-U0)  t  . 


(B-36) 


This  process 


V0~U0 


continues  until  the  length  L  is  reduced  to 
The  maximum  penetration  occurs  at  t„»x  and  is: 

p  ■  r 

V0-U0  0  • 


zero  at  a  time 


(3-36) 


At  any  instant  of  time  while  penetration  occurs,  we  have: 

y°~u°p(t)  +l(  t)  =  l0  . 


( B-37 ) 


The  exact  solution  (solved  for  in  the  manner  of  Equation  B-10)  and  the 
approximate  solution  (refer  to  Equation  B-12)  are: 


vo~Uo  _  *o+\/V2+<*o(l-V2)  m  H  + 

U0  l-a0  l-a0 


where  a0  =  H/(kp  pp  V02)  .  E<juation  B-38,  with  kp  =  1/2,  leads  to  the  Si  equation 
(Equation  16)  of  Section  III  which  is: 


SxiZa) 


foz  ZC<Z0£°° 


(B-39) 


and  S1(Z0)=0  for  0<Zo<Zc=n  . 

Some  additional  insight  into  the  penetration  process  can  be  obtained  without 
integrating  the  system  of  differential  equations.  From  Equation  B-23: 

dP=-ZR.dL  ,  ( B-40 ) 

independent  measurements  of  AP,  AL,  U  and  V  can  be  used  to  deduce  the  values  of 
H  -  Y  from  experimental  observations  assuming  that  the  penetration  process 
conforms  approximately  to  Equations  B-l  through  B-4.  Note  that  the  H  -  Y  value 
so  determined,  depends  also  on  the  shape  factors  kp  and  kt.  The  values  for  the 
shape  factors  can  sometimes  also  be  deduced  from  experimental  observations. 

Another  useful  relation  can  be  deduced  for  the  crater  diameter  (ignored  in 
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this  one  dimensional  treatment)  by  invoking  the  work-energy  hypothesis.  To  an 
approximation,  the  energy  lost  by  the  rod  is  partitioned  into  work  on  the  target 
(i.e.,  making  a  hole  of  diameter  Dt/D  =  h  >  1)  and  work  to  "fluidize"  the  rod. 
The  energy  balance  at  any  instant  of  time  is: 

—^(±ppApLV2}  =  AcHu+ApY  (V~U)  .  ( B-41 ) 

With  At/Ap  =  h2  >  1  and  using  Equations  B-2,  B-3,  we  obtain: 

YV*-^p[HV-U)V2^h2HU+Y(V-U)  .  ( B-42  ) 

Using  Equations  B-8  and  B-12,  this  can  be  expressed  as: 

h2  =  l+-^/-t^-2k-a\  .  ( B-43 ) 

2 H  \ 1-a  p  I 

For  the  special  case  of  kp  =  1/2,  this  becomes: 

h2  =  14  A  PPV2f  .  ( B-44 ) 

2  H  \  1-a  ) 

which  reduces  to  the  well  established  experimental  relation  of  h  being 
approximately  linearly  dependent  of  V  at  very  high  velocities  (a  «  1). 


What  has  been  assumed  so  far  is  that  H  and  Y  are  constant.  Some  materials 
undergo  large  deformations  in  a  rate  dependent  mode  -  they  exhibit  visco-plastic 
behavior.  The  H  and  Y  values  also  depend  on  the  shape  of  the  penetrator/target 
interaction  region  that  is  only  approximately  described  by  the  simple  shape 
factors  kp  and  kt.  That  Y  can  be  treated  as  a  constant  is  not  obvious.  As  long 
as  the  crater  cross-section  made  in  the  target  is  sufficiently  large  to  allow  for 
radial  "flow"  of  material  away  from  the  erouiag  pc**e^raLoj. ,  une  Y  value  is 
expected  to  be  about  equal  to  the  stress  of  a  penetrator  cylinder  loaded  on  an 
unconfined  end,  resulting  in  large  deformations  on  that  end  (radial  flow,  "mush¬ 
room"  formation).  An  estimate  of  the  value  of  H  can  be  determined  from  cavitation 
analysis  (e.g.,  Hanagud,  Ross'  or  Hopkins”). 

Despite  the  limitations  imposed  by  the  assumptions,  this  simple  LRP  model 
is  quite  useful. 
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APPENDIX  C 

EXAMPLES  OP  PARAMETRIC  VARIATION 
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The  figures  in  this  appen  ’ ■' x  are  generated  by  using  the  following  equations. 

P(Z0. 0)  =  (L-mD)  Si(Zo,0)+mD  SD(Z0,0)  ,  (C-l) 

S3(Z0 ,  0)  =  -jj|l - + -  *—  —  j  for  p<a0<~  ,  (C-2) 


o 

II 

o' 

o 

for  0<a0£n  , 

r2l  2  +  p  \ 

2 Y  \  1  +  p/ 

H 


2  +  p 


SD(Z0 ,  Z)  = 


Z=  0 


and 


£j. 

Pp 


(C-4) 


(C-5) 


(C-6) 


(C-7) 


In  these  equations,  P(ZQ,0)  is  the  depth  of  penetration  into  the  target,  L  is  the 
initial  length  of  the  penetrator,  m  is  an  empirical  constant,  D  is  the  penetrator 
diameter,  pt  is  the  target  density,  pp  is  the  penetrator  density,  H  is  target 
material  resistance  to  penetration  (units  of  pressure)  and  VQ  is  the  impact 
speed. 

In  the  following  figures,  the  target  density,  p{,  is  7850  kg/m3,  H  =  5  GPa 
and  m=2/3.  The  constants  kt,  kp  and  Jcq  of  Appendix  A  and  Appendix  B  are  each  set 
to  1/2.  The  L/D  ratio  (denoted  by  f)  takes  on  the  vaules  of  1,  2,  5,  10,  20,  50 
and  100  for  each  respective  curve.  For  each  parameter  plotted  on  the  ordinate 
there  are  a  pair  of  figures.  The  first  of  the  pair  is  for  pp  =  17300  kg/m3 
(corresponding  to  91%  tungsten  alloy)  with  fi  *  0.6736  and  the  second  is  for  pp 
=  18600  kg/m3  (corresponding  to  DU-3/4Ti)  with  p  =  0.6496.  The  abscissa  in  all 
figures  is  impact  velocity.  The  ordinate  for  each  pair  is  as  follows: 
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Figures  C-l  and  C-2:  Penetration  per  unit  length  (P/L)  evaluated  from 

|  =  -g^(Z0)+(1-^)s3(Z0)  , 

Figures  C-3  and  C-4:  Penetration  per  unit  diameter  (P/D)  evaluated  from 

^=mSD(Z0)+  (f-m)S3(Z0)  , 

Figures  C-5  and  C-6:  Penetration  per  cube  root  of  unit  volume  (P/VOL1/3) 

evaluated  from 

— =(-£)"* js3(Z0)l  , 

VOL3 


Figures  C-7  and  C-8: 


Penetration  per  cube  root  of  unit  kinetic  energy  (P/E1/3) 
evaluated  from 


and 

Figures  c-9  and  C-10:  Penetration  per  cube  root  of  unit  mass  (P/M1/3)  evaluated 

from 
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FIGURE  C-2 


penetration  /  DIAMBTER 


C  -  6 


C  -  7 


No  of 

No  Of 

CoDies  Oraanization 

CoDies 

Oraanization 

2 

Administrator 

1 

Commander 

Defense  Technical  Info  Center 

US  Army  Missile  Command 

ATTN:  DTIC-DDA 

ATTN:  AMSMI-RD-CS-R  (DOC) 

Cameron  Station 

Redstone  Arsenal,  AL  35898-5010 

Alexandria,  VA  22304-6145 

1 

Commander 

1 

HQDA  (SaHu-iR) 

US  Army  Tank-Automotive  Command 

WASH  DC  20310-0001 

ATTN:  AMSTA-TSL  (Technical  Library) 

Warren,  Ml  48397-5000 

1 

Commander 

US  Army  Materiel  Command 

1 

Director 

ATTN:  AMCDRA-ST 

US  Army  TRADOC  Analysis  Command 

5001  Eisenhower  Avenue 

ATTN:  ATRC-WSR 

Alexandria,  VA  22333-0001 

White  Sands  Missile  Range,  NM  88002-5502 

1 

Commander 

(Class,  only)  ■) 

Commandant 

US  Army  Laboratory  Command 

US  Army  Infantry  School 

ATTN:  AMSLC-DL 

ATTN:  ATSH-CD  (Security  Mgr.) 

Adelphi,  MD  20783-1145 

Fort  Benning,  GA  31905-5660 

2 

Commander 

(Unclass,  only)  ■] 

Commandant 

US  Army,  ARDEC 

US  Army  Infantry  School 

ATTN:  SMCAR-IMI-I 

ATTN:  ATSH-CD-CSO-OR 

Picatinny  Arsenal,  NJ  07806-5000 

Fort  Benning,  GA  31905-5660 

2 

Commander 

1 

Air  Force  Armament  Laboratory 

US  Army,  ARDEC 

ATTN:  AFATL/DLODL 

ATTN:  SMCAR-TDC 

Eglin  AFB,  FL  32542-5000 

Picatinny  Arsenal,  NJ  07806-5000 

Aberdeen  Provina  Ground 

1 

Director 

Benet  Weapons  Laboratory 

2 

Dir,  USAMSAA 

US  Army.  ARDEC 

ATTN:  AMXSY-D 

ATTN:  SMCAR-CCB-TL 

AMXSY-MP,  H.  Cohen 

Watervliet,  NY  12189-4050 

1 

Cdr,  USATECOM 

ATTN:  AMSTE-TD 

1 

Commander 

3 

Cdr,  CRDEC,  AMCCOM 

US  Army  Armament,  Munitions 

ATTN:  SMCCR-RSP-A 

and  Chemical  Command 

SMCCR-MU 

ATTN:  SMCAR-ESP-L 

SMCCR-MSI 

Rock  Island,  IL  61299-5000 

1 

Dir,  VLAMO 

ATTN:  AMSLC-VL-D 


1  Director 

US  Army  Aviation  Research 
and  Technology  Activity 
ATTN:  SAVRT-R  (Library) 
M/S  219-3 

Ames  Research  Center 
Moffett  Field,  CA  94035-1000 


23 


Intentionally  left  blank. 


24 
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for  which  the  report  will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source 
of  ideas,  etc.) _ _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars 
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